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ABSTRACT 


The magnetic field measurements made by the magnetic field 
experiment on the IMP 3 (Explorer 28) spacecraft have been examined 
at the time of geomagnetic s^s.c. events. 36 such events occurred 
while IMP- 3 was in the interplanetary medium during 1965, 66 and 67 
and have been analysed. Of these events 8 must have been tangential 
discontinuities, 2 are either tangential discontinuities or 
rotational discontinuities and 26 are possible shock waves. 

(2 of these 26 events have been shown by other authors to be 
shocks. ) These 26 possible shocks have similar magnetic 
signatures; an increase of 20 $ or more in the magnetic field 
magnitude and a relatively small (always less than 90°) change 
in direction. The larger s.s.c. events were more likely to be 
caused by possible shocks while the smaller events were often assoc- 
iated with tangential discontinuities. The orientation of the dis- 
continuity surfaces of the 26 possible shocks shows a preference to be 
aligned somewhere between a direction perpendicular to the sun-earth 
line and a direction tangent to the local spiral angle of the 
magnetic field. It was possible to associate solar flares with 
lk of the 26 possible shock events. Of these Ik a reliable 
orientation was deduced for 8 events. By considering the 
orientation of these 8 events in relation to the position of 
the parent flares on the solar disk it is suggested that a. typical 
shock front propagating out from the sun has a radius of 
curvature less than but of the order of 1 a.u. 



INTRODUCTION 


A number of authors have considered discontinuities in 
the solar wind associated with geomagnetic disturbances 
(Gold, 1955 } Sonett et al., 1964 ; Hirshberg, 1965 J Col- 
burn and Sonett 1966 j Gosling et al., 1967a, ' 1967b, 

1968 ; Van Allen and Ness, 1967 ; Ogilvie et al., 1968 ). 

It is agreed that geomagnetic sudden commencement (s.s.c.) 
events are caused by a discontinuous increase in the plasma 
pressure on the magnetosphere. Satellite studies of these 
magnetohydrodynamic discontinuities have revealed that some 
are shockwaves (Sonett et al. , 1 964 ; Gosling et al., 1967a; 

1968 , ) and some are tangential discontinuities (Gosling et al. # 

1967b ). 

These s.s.c. 'associated discontinuities are thought to be 
the boundary between plasma ejected from the sun in a solar 
flare and the quiet, pre-flare solar wind. The propagation of 
the flare plasma out from the sun has been considered by 
various authors (Gold, 1959 ; Parker, 1963 ; Hirshberg, 1968 ) 

and various models have been suggested. The measurement of the 
magnetic field associated with this type of s.s.c. related event 
gives information about the local orientation of the discontinuities 
and thus can help to determine the large scale structure of these 
features. In this paper a number of sudden commencement-associated 
events have been studied with the magnetometer on the IMP 3 (Explorer 28) 
spacecraft. Two findings of this study can be summarized as follows. 


1. Both shock waves and tangential discontinuities are 
capable of causing an s.s.c. event. The larger 
events are more likely caused by shock waves. 

2. When it passes the earth a typical shock front caused by 
the flare- associated outburst appears to have a radius 
of curvature of less than, but of the order of, 1 a.u. 

The measurements from 3 6 different s.s. c.- associated events which 
occurred while the IMP 3 spacecraft was in the interplanetary 
medium between June, 1965 end January, 1967 will be presented and 


discussed. 


THE EXPERIMENT 


The IMF 3 spacecraft was launched on May 29, 1965 into a 
highly eliptical orbit. Initially its apogee was hi. 9 Re 
(Re*6378 km) and decreased to 36.3 Re during the spacecraft's 
lifetime of approximately two years. The perigee varied from 
only slightly more than 1 Re (a few hundred km altitude) at 
launch to as high as 6. 5 Re late in its lifetime. The initial 
sun- earth-apogee angle was 121 degrees. The orbital period was 
initially 140 . 6 hours or 5.86 days .and remained within 1.5 
houis of this value during the lifetime of the satellite. 

During the months of June through January of each year the 
spacecraft spent at least some time in the interplanetary 
medium, upstream from the bow shock. During the months of 
February, March, April and May IMP 3 never crossed the bow 
shock and therefore did not sample the interplanetary medium 
directly. 

The magnetic field experiment on IMP 3 was identical to 
that used on IMPs 1 and 2. It consisted of two monoaxial 
spinning fluxgate magnetometers. The instrument has been 
described previously in detail by Ness et al. (196h) and by 
Fairfield and Ness (1967). 

Unfortunately one of the magnetometers failed during the 
launch sequence boom erection so that measurements are obtained 
only half as frequently, that is every kl seconds instead of 
every 20.5 seconds as on the earlier IMP 1 and 2 spacecraft. 


There is a gap of 123 seconds after every 6 data points when the 
rubidium magnetometer is read out. Each measurement is made in 
a period of 4.8 seconds by making • least squares fit of the 
spin modulated output of the detector to a sine wave. The root 
mean square (RMS) deviation of this fit is used as an indication 
of the goodness of the fit. Detailed analyses of the possible 
errors in this scheme have been published by Fredricks et al. ( 196 /) 
end Fairfield and Ness ( 1967 ). 

Direct comparisons with other spacecraft (Pioneer 6 and 7 
and Explorer 33) throughout the IMP 3 lifetime give excellent 
agreement. Such comparisons indicate that the absolute error of 
the instrument is probably l)f (C»10“5 gauss) or less. The 
relative accuracy important in this study of discontinuous changes 
in the field, is even better, especially over the short time 
periods of interest in this work (of the order of minutes). 


ANALYSIS 


Events were chosen for study by consulting the preliminary 
report of sudden commencements published quarterly in the 
Journal of Geophysical Research prior to 19 66 and subsequently in 
Solar- Geophysical Data* at the same 3 month intervals. Geo- 
magnetic impulses identified by 10 or more magnetic observatories 
as s.s.c. events ("Sudden commencements followed by a magnetic 
storm or a period of storminess" ) were selected for months of 
June, 1965 to January, 1966 and June, 1966 to January, 1967. 

Then the IMP 3 magnetometer records were examined near the 
time of each reported s.s.c. For the events which occurred 
while IMF 3 was in the interplanetary medium the associated 
discontinuity in the magnetic field was selected and a computer 
program used to calculate the average field on each side of the 
discontinuity. Periods of between 5 and 15 minutes were chosen 
over which to average, the interval length being chosen so that 
the field was reasonably steady throughout. This meant that the 
averages consisted of between 6 and 20 individual measurements. 

In addition to the field averages B^ and and the 

associated RMS deviations, various quantities were calculated 


including : 

the vector change, 

the shock normal, n, =(? x §) xa! / \ c?x ■? ) xifel 
the normal component, B*. « ? • i • & 

the parallel components, ** - II B* , and - * % ; etc. 


* Published monthly by the Environmental Science Services 
Administration of the U. S. Department of Commerce. 



From these quantities it was determined whether or not 
the discontinuity might be a shock. The conditions for a fast 
shock (discussed by Colburn and Sonett, 1966 ) are the following: 

B 6 > ‘ B l T 2 >Tl and * ^>, 0 . Because plasma data are 

not available on IMP 3 (the plasma experiment failed) the 
conditions on density, , and temperature, T, cannot be checked, 
but the requirement on the fields was checked. The events with no 
clear increase in the field magnitude, B, or with Bi • "B* 4 0 were 
identified as tangential discontinuities (or if B 1 = B* possibly 
a rotational discontinuity) and the other events as possible 
shocks . 

The calculation of Zb and of n are used to give information on 
the orientation of the surface of the discontinuity. Ideally n is 
all that is needed, but if the discontinuity is not a shock, then 
n as calculated does not give the normal direction. In all 
cases, however, Zb must lie in the plane of the discontinuity 
since the flux across the surface of the discontinuity is conserved. 
In addition the calculation of n becomes inaccurate when B L is 
nearly parallel to . Also the calculation of n is dependent 

on the instrumental zero levels while Z? is not. For these reasons 
it seemed useful to use the vector AB as an indication of 
orientation of the discontinuity surface rather than the unit 
normal A. If the normal to the surface is assumed to lie in the 
ecliptic plane then a knowledge of is sufficient to determine 
the orientation uniquely. 


Another advantage of using /il is that the RMS error in is 


easily calculated. It can be shown to be given by the Pythagorean 
sum of the RMS errors in the individual components of 1?- and 3^ . 
Thus it is a simple matter to estimate the accuracy of the measure- 
ment of & by considering the quantity where O' is this RMS 

error in Zt. If this quantity is greater than 1 then the measure- 
ment is rather poor and the uncertainty in the direction of AB 
is large. On the other hand, if this quantity is less than 0.5 
then the directional error in is statistically expected to be 
less than 1 3 <f . 

Another characteristic of the discontinuity, its velocity 
in the vicinity of the earth, can be calculated from the observations 
by dividing the perpendicular distance between the discontinuity 
surface and the earth by the delay time between the recording of 
the event at the spacecraft and its record at the earth. Unfortu- 
nately because of the different propagation characteristics of the 
magnetosheath and the magnetosphere the average velocity obt ain ed 
in this way is not necessarily representative of the surface in the 
interplanetary medium. Alternatively this uncertainty may be 
viewed as an uncertainty in the time at which the surface passes the 
earth. A recent estimate of the propagation time within the 
magnetosphere yields values of the order of 90 seconds (Sugiura, 

1965 ). The propagation time through the magnetosheath 

depends on the direction of arrival and is estimated to be of 
the same order as o r slightly greater than this time of 



90 seconds. Thus there is an uncertainty of 3 minutes or more 
in the time of arrival of the surface at the earth. For delay 
times that are long compared to 3 minutes this uncertainty 
vould he unimportant, hut for short delay times it introduces 
large inaccuracies. Because of this velocities were 
only calculated for those shocks with delay times greater 
than 6 minutes. 

Finally, the relation of the s.s.c. events to solar flares 
was investigated. All flares of. importance It or greater 
occurring between 10 and 120 hours before the s.s.c. were 
considered and the largest disturbance (or the flare with 
associated Type II radio emission) was picked as the most likely 
parent flare. If there were several flares of comparable 
magnitude or no flares greater than importance 1 then no 


identification was made. 


THE OBSERVATIONS 


In all there were 52 s.s.c. events identified by 10 or 
more magnetic observatories during the months of June, 1965 
to January, 19 66 end June, 1966 to January, 1967# 0 f these,, 36 

occurred while IMP 3 was in the interplanetary medium. The 
magnetometer observations from IMP 3 for these 36 events are 
summarized in Table 1. 

The first column of this table is simply a sequential 
numbering of the events. They are listed in chronological 
order. The next three columns show the date, universal time, 
and the number of magnetic observing stations reporting the 
evenly, respectively. The number of stations reporting the 
event is used as a statistical, estimate of the significance 
of the event. Its usefulness is illustrated by the observation 
that each of the 9 events reported by more than 50 
stations was found to be a possible shock, while of the 9 
events reported by less than 20 stations only 4 were possible 
shocks . 

The rest of the columns in the table relate to the IMP 3 
measurements. The column headed ORB is the IMP 3 orbit number. 
When two or more events occurred during the same orbit the 
events are distinguished by a letter (e.g.,the two events on 
orbit 7 are labelled 7a and 7b). These orbit number 
identifications will be used on subsequent figures for identi- 
fication purposes. The column headed UTiat is the time of the 
passage of the discontinuity at IMP 3 along with the uncertainty 



in this time. In all but one case the uncertainty is just half 
the time between successive measurements -™that is, 20.5 
seconds (about 0.4 minutes) if the jump occurred between 2 
successive points 4l seconds apart or 6 1. 5 seconds (about 
1.0 minutes) when the jump occurred across the gap where the 
rubidium magnetometer is read out. In the one unusual case 
(event 9) the field underwent first a decrease and then 5 
minutes later a sharp increase and it was unclear just what 
time should be selected, hence the large At. 

The next three columns give the cartesian components of 
the spacecraft position vector in geocentric solar ecliptic 
coordinates. Xse is toward the sun, Zse toward the north 
ecliptic pole and Yse in the dusk meridian forming a right- 
handed system. The next 6 columns list the average magnetic 
field components along with the RMS deviation of the average. 

"3^ is the field before the discontinuity and is after. 

Again the coordinates are solar ecliptic. The next three 
columns give the components of the vector change in the field, 

AB. The column labelled °7 ab is the ratio of the RMS deviation 
of the vector change to the magnitude of the vector change and 
measures the reliability of the measurement. As mentioned in 
the preceding section, if this number is 0. 5 or less the 
measurement is good. The last column, labelled Sjf, is a 
sequential numbering of the possible shocks. The events for 


which 4 B 1 or * Bf4 0 have no entry in this column. 

Thus it is seen that of the 36 s.s.c. events which occurred while 
IMP- 3 was in the interplanetary medium, 26 of them were possible 
shocks. 

Examples of the two types of events are shown in Figure 1. 
Figure 1a, is event number 1 in Table 1. It occurred on June 8, 

1965 and is a possible shock. The discontinuous increase in the 
field magnitude and the relatively small change in field 
direction are characteristic. The event shown in Figure lb is 
a tangential discontinuity. It is characterized by a large 
change in the field direction and thus cannot be a shock. This 
is event 3 in the table and occurred July 6, 1965. 

In Figures 2 and 3 information from all 26 possible shock 
events has been superimposed for the purpose of comparison. 

Figure 2 shows the 26 possible shock events projected onto the 
ecliptic plane (by rotation about the earth-sun line). The 
line plotted is the projection of on the ecliptic and 
represents the orientation of the event if it is assumed that 
the discontinuity surface is normal to the ecliptic or equivalently 
that the surface normal lies in the ecliptic. The error in each 
measurement may be obtained from the °/ab column of Table 1. 

The arrow plotted for each event is the projection of the unit 
normal on the ecliptic. This normal is calculated assuming 
that the discontinuity is a shock. This figure shows that the 


orientation of the discontinuity surfaces is not in general 
perpendicular to the sun-earth line as is often assumed. 

Figure 3 shows that in en average sense there is some 
order to the orientation of the discontinuities. By plotting 
the direction of the change vectors in this way it is 
evident that there is a preferred orientation for these vectors. 
There is apparently no preference in the 9 direction, hut in the 
azimuthal angl$ If , the events scatter about an angle of about 
115° (or 295°). It is probably significant that this direction 
is intermediate between the perpendicular to the sun- earth line 
(90° or 270°) and the magnetic field spiral angle (135° or 315°). 
The dashed lines are the average angles for the outward and 
inward directed groups of changes. 

Thus statistically there is evidence for a preferred 
orientation of the surfaces of the discontinuities responsible 
for s.s.c. events on the ground. Furtner information can be 
obtained by examining the events for which parent flares could 
be tentatively selected. 

The 26 possible shock events together with information 
about the identification of the probable parent flare are 
listed in Table 2. The first two columns are the date and time of 
the event on the ground as in Table 1 and the third column 


gives the IMP 3 orbit number. The next 5 columns give for the 
flare: the date, the time of onset, importance, latitude 
and longitude^ respectively. Note that for the first event 
there was a type II radio burst but no flare observed. 

The next two columns give the delay time between the 
occurrence of the flare and the onset of the s.s.c. 
at earth and the average sun-earth transit velocity implied Tpy 
this delay. The last columns gives the estimate of the 
near earth velocity calculated from the IMF 3 measure- 
ments for which the time delay was longer than 6 
minutes . There is an estimated uncertainty of 50$ associated 
with this last velocity. 

Tills table shows that flares were identified for lU of 
the 2 6 possible shock events. It is interesting that the 
identification of flares for only about half of the 2 6 possible 
shock events is statistically similar to the findings of 
Kuleshova (1961) who found geomagnetic activity associated 
with 50$ of the class three flares which 



occurred over a three year period. Of the remaining 12 events 
for which no flare could be found 9 had no flares greater than 
importance 1 in the period, between 10 and 120 hours before the 
s.s.c. while the other 3 had multiple flares of comparable 
importance. Examination of the IMP 3 measurements related to 
the ll flare- associated discontinuities showed that for 6 of 
these events the RMS deviations were less than 65$ of the 
magnitude of Al. Two others were included because plasma data have 
already been published for these events (January 20 , 19 66- 
Gosling et al., 1968; and July 8, 1966 —Lazarus and Binsacl^ 
i960 ). For the January 20 event 0 "jlB was 0.8 and for the 

July 6 event additional data from the Explorer 33 experiment 
(Ness and Taylor, 1968) was used to supplement the IMP 3 
measurements. Thus reasonable orientation data is available 
for a total of 8 flare- associated, events. These are plotted 
in F igure 4 , in an ecliptic plane view of the earth-sun system. 

The discontinuity vector, projected on the ecliptic is 

plotted at the heliocentric longitude of the earth at the time 
of the flare. The longitudes are specified as seen from the 
earth. For example the flare associated with the s.s.c. of 
July 15, 1966 during orbit 71 occurred at 90 ° W on the solar 
disk and thus is plotted at 90 ° W in the figure. This method 
of plotting the events strongly suggests that for the events 
occurring west of the central meridian on the sun the IP angle 


of the discontinuity surface is less than 90° while those 
occurring East of the central meridian have a (f angle greater 
than 90 °. This is in agreement with the work of Hirshberg 
(1968) and is also predicted by the qualitative picture of 
the plasma outburst from a flare presented by Gold ( 1959 ). 

A final observation can be made about the velocity of 
the shock in the vicinity of the earth relative to the 
average velocity of propagation between the sun and the earth. 
Although it was possible to determine the near earth velocity 
with reasonable accuracy in only 2 of the 14 flare-associated 
events, both cases indicate that the near earth velocity is 
less than the average sun to earth velocity in agreement with 
the findings of Gosling et al. (1968) and of Ness and Taylor 
(1968). 


DISCUSSION 


The identification of the events as to the type of 
discontinuity is not unambiguous from the magnetic field data 
alone. It is possible to rule out shocks in cases where the 
directional change is too large. It Is also possible to rule 
out rotational discontinuities whenever the magnitude of the 
field changes. However, tangential discontinuities overlap 
both shocks and rotational discontinuities in their possible 
magnetic signatures end thus cannot be ruled out without plasma 
data. In the 36 events studied in the work reported here 10 
were clearly not shocks. Of these 10, two might possibly have 
been rotational discontinuities but the rest must have been 
tangential discontinuities because the magnitude of the field 
changed by more than 20$. Among the 26 possible shocks were 2 
events for which plasma data have been published (January 29, 
1966 - Gosling et al., I96Q' and July 8, 1966 -Lazarus and Binsack^ 
1968 ). These two events were clearly shocks and thus it 

does not seem unreasonable to assume that many of the remaining 
2b possible shock events were indeed shocks. As a rough 
estimate, perhaps half of the 3 6 events studied were 
shocks and the other half tangential discontinuities. 

Of the 26 events labelled as possible shocks a probable 
flare can be associated with only lh of them' others were 
ambiguous. The events which were related to their parent flares 
showed orientations in good agreement with the findings of 



Hirshberg (1968). A simplified picture of her resulting 
large scale structure for the shock is included in Figure 4 
as the dashed line. This line is an arc of a circle of radius 
0.75 a.u. and centered on the 0° 1 ine 0.5 a. u. from the sun. It 
vas drawn as an aid in visualizing the large scale structure 
of the interplanetary shock. 

It seems clear that the measurements imply that the typical 
shock front has a radius of curvature less than, but of the order 
of, 1 a.u. Such a picture seems to be intermediate between 

Parker’s (1963) spherical blast wave model and Gold's (1959) magnetic 
tongue model. It is not really inconsisle nt with either. 

A close look at figure 4 reveals two events which are not 
consistent with the shock surface drawn. These two events can 
be understood if one is viewed as an example of Gold's model 
while the other is an example of Parker's. Orbit 101a shows 
a discontinuity vector with a large component in the sunward 
direction. If this discontinuity were the result of a magnetic 
tongue reaching out from the central meridian, perhaps deflected 
slightly along the spiral angle^ then the orientation measured 
would be expected. Orbit 101b is also at variance with the shock 
front drawn because it has an azimuthal angle of slightly more 
than 90° where it should be less. It is, however, not greatly 
different from 90° end thus might be interpreted as an example 
of a spherical blast wave. Of course for flares near the 
central meridian the three pictures (Gold's, Parker's end 
Hirshberg's) all predict the same orientation. 


Because the events on the Eastern limb of the sun seem to 
show a slightly larger average inclination while apart from 
orbit 101a those on the Western limb show slightly less 
inclination than the shock drawn, there seem to be a suggestion 
that the best center of symmetry might be to the west of the 
central meridian. This shift would be slight (*>10 o ) and the 
data cannot be said to be conclusive. Hirshberg's study 
seemed to show a similar skewness. 

Another conclusion with marginal support in the data is 
related to the relative orientation of the shock and tangential 
discontinuity events. Specifically, the shock events might be 
expected to scatter about an orientation of 90° while the 
tangential discontinuities might be found to scatter about the 
spiral angle direction or about 135°. Then with roughly equal 
numbers of each type with similar scatter about their respective 
means, the resulting distribution would be very much like that 
shown in figure 3« The average of 'the change vector would 

be between 90° and 135° (or 270° and 315°) end few or no points 
would lie in the rang^k 0° - 45® and 180° - 225° precisely as 
the figure shows. Clearly this is highly speculative but it 
does seem to be a plausable explanation of the data. 



CONCLUSIONS 


In this study 3 6 sudden commencement events which occurred 
while IMP 3 was in the interplanetary medium have been examined. 

It has been found that there are two distinct types of 
magnetohydrodynamic discontinuity associated with these events. 

Some are tangential discontinuities while others are shock 
waves. It is suggested that among the s.s.c. events reported 
by more than 10 stations both types may be found. The events 
identified by a large number of stations are more likely to 
be shocks than those reported by relatively few stations. 

The possible shock events are characterized by a discon- 
tinuous increase in the magnetic field magnitude (20 $ or more) 
accompanied by a small change in the field direction (always 
less than 90°). The events cannot positively be identified as 
shocks without information about the density and temperature of 
the plasma. Some events can, by elimination, be positively identi- 
fied as tangential discontinuities. They are characterized by a 
relatively large directional change and a significant change in 
the field magnitude. 

When all the possible shock events are taken together a 
preferred azimufchal orientation is found between 90 ° and. 135 ° 

(or 2T0° and 315°). When the 8 reliable measurements of events 
for which the parent flare is identified are plotted as a 
function of heliocentric longitude it is apparent that a shock 
front with a radius of curvature of the order of, but less than 


1 a.u. is suggested in agreement with the findings of Hirshberg 
(1968). There are indications that some events may be of 
Parker ' s blast wave type while others may be like Gold ' s 
magnetic tongue model, but apparently the bulk of the events 
are intermediate between these two models. 

Although useful information can be obtained from the 
magnetic field data alone it is clear that the complete picture 
cannot be obtained without complete plasma data. In addition, 
it would be useful to have simultaneous and independent data 
from other spacecraft in cis-lunar space. An ideal opportunity 
for such a study exists during 1967 and 1968 when the three 
spacecraft Explorers 33, 3 U and 35 are all operating and much 
of the time they are making simultaneous measurements of both 
the magnetic field and the plasma at three separate points 
in space. Studies of interplanetary discontinuities using 
simultaneous data from all three spacecraft should provide 
definitive information on these interesting events. 
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FIGURE CAPTIONS 


Figure 1 . 


Figure 2 . 


Figure 3. 


Figure 4 . 


The magnetometer measurements from IMP - 3 vhile 
in the interplanetary medium along with ground 
magnetograms. The field magnitude scale is in 
gammas. Ais a measure of the short period 
( ~ 5 sec. ) fluctuations in the field. 0 is the 

latitude angle of the field (90° is toward the north 
ecliptic pole, 0° is in the ecliptic plane) and is 
the azimuth angle. (0° is toward the sun, 90° 
toward dusk etc.) The coordinates are geocentric 
solar ecliptic. 

A plot of the orientation of the observed discontinuity 
surfaces relative to the earth and its magnetosphere. 

Each surface is labelled with the IMP - 3 orbit number for 
identification purposes. All events actually occurred in 
the interplanetary medium although some are inside the 
average position of the bow shock. 

A scatter plot of the direction of AB. As in Figure 1 9 
is the latitude angle of the field and ty is the 
azimuth angle. 

A plot of the orientation of 8 probable shock surfaces 
at the appropriate heliocentric longitude relative to the 
flare. The dashed line is an arc of a circle of radius 
0.75 a.u. centered on the 0° line 0.5 a.u. from the sun. 
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